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(b) 

Figure 1. Physical scatterers of the same scattering length but different radii 
introduced into the rectangular billiard. Black circle shows the scatterer, the radius of 
the dashed circle is equal to the scattering length. 



(a) 



1. Introduction 

The singular perturbed square billiard, also called Seba billiard [lj by several authors, 
is one of the key models for quantum chaotic systems. While the unperturbed square 
billiard with side ratio l/(\/5 — 1) shows Poissonian level-spacings statistics [2] and the 
Sinai billiard, being proved to be a fully chaotic system [3], exhibits GOE-statistics [I], 
the "intermediate" case of a singular perturbed billiard is expected to demonstrate some 
transient behavior jH [51 [7J . However it was reported, that the billiard with a point 
perturbation can exhibit "fully developed quantum chaos" [HE]. It may seem strange, 
since a point perturbation has almost no influence on the classical phase space of the 
billiard. Thus it is natural to assume that in the semiclassical limit the billiard with a 
point perturbation shows a similar statistics to an unperturbed one. 

The proposed explanation of the given paradox was based on the argument that for 
zero-range perturbation for any wavelength one can never reach the limit of the classical 
billiard with a point-perturbation, since the wavelength is finite while the perturbation 
radius is zero. Therefore the quantum system becomes chaotic while its classical analog 
is almost integrable. This argument seemingly was justified experimentally [6]. 

In the presented report we show that the level-spacings statistics of Seba billiards 
actually tends to Poissonian when the number of taken eigenvalues tends to infinity. 
These findings are in accordance with the intuitive "classical" argument given above, 
but are in contradiction with previous theoretical [U |7j and experimental [6] results. For 
a narrow window of eigenvalues some conclusions of Refs. [U [7J remain valid, however 
one cannot directly apply the results to a wide eigenvalue range. This discrepancy 
traces back to the procedure applied by the authors to take care of the singularity of 
the Green function, by replacing the "bare" coupling constant by a renormalized one 
absorbing the infinity. Although this renormalization technique is standard in quantum 
electrodynamics, it is not appropriate to compute the spectrum of a Seba billiard. In 
this paper we present the suitable renormalization procedure. 

Let us turn to the experimental microwave Seba billiard [6] and explain why the 
interpretation of the obtained distributions was erroneous. Theoretically in all cases the 
scatterer was treated as a point scatterer, which means physically that the characteristic 
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Figure 2. Sketch of the level-spacings statistics evolution computed for a given 
fixed number of eigenvalues of a billiard with a small perturbation whose radius a is 
significantly smaller than its scattering wavelength j3. The direction from left to right 
corresponds to increasing eigenvalues. The "semi-Poissonian" statistics (a) corresponds 
to the range a <C j3 < A, the "Poissonian" statistics (b) corresponds to the range 
a <C A <C (3 and the "GOE" statistics (c) corresponds to the range a ~ A <C (3. Here A 
is the characteristic wavelength. Here A is the characteristic wavelength. The names 
of distributions are written in parenthesis since the corresponding plotted curves keep 
the essentials of these statistics but may differ from real distributions. 

wavelength A of the field inside the cavity is much larger than the radius a of the 
scatterer. At the same time, the scattering length (3 of the given scatterer is, generally 
speaking, a free parameter depending on the internal structure of the scatterer (e. g. 
given material of the coating, radius of the metallic core etc). We show below that the 
influence of the point scatterer is significant when A > (3 and vanishes when A <C (3 in 
accord with a classical limit. To treat the scatterer as a point perturbation we obligatory 
should require o < A. Combining the last two estimations we conclude that to cover 
experimentally the classical limit of a Seba billiard one needs to create a scatterer with 
a <C (3 (see figure [□ (b)). 

However, experimentally one often has the situation shown in figured] (a), i.e. 
a ~ (3. This means that in the regime A ~ (3 the corresponding billiard should be 
treated as a quantum Sinai billiard, but not a Seba billiard. In this case the "classical" 
limit with a point perturbation can not be achieved. 

The expected "experimental" evolution of the level-spacings statistics computed 
for a given number of resonances taken at different frequencies is plotted in Fig. [2j 
It has been assumed that the radius of a small scatterer a <C (3, but remains finite. 
Figures 2(a), 2(b) correspond to the Seba billiard approximation when the radius of 
the perturbation can be neglected. Figure 2 (c) shows the "GOE" statistics of the Sinai 
billiard in the regime where the wavelength of the field is comparable with the radius 
of the perturbation. 

In what follows we take the limit a — * which means that figure 2 (c) can not be 
reproduced within the framework of the considered approach. Thus we restrict ourself 
to the mathematical model of a point perturbation as it was done by the authors of Refs. 

mm- 
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2. Point perturbation of the billiard 



Let us now come to the theoretical and numerical study of the Seba billiard and show 
why the previous treatments have given incorrect results. Following [HI H] let us first 
introduce a point perturbation of the billiard at point R. Basically we will construct the 
"self-adjoint extension" of the unperturbed "Hamiltonian" . This approach was already 
previously used in [10J. 

For the unperturbed billiard the eigenfunctions ip and eigenvalues k 2 = k\ obey the 
equation 

(A + k 2 )ip(r) = 0. (1) 

If ip(R) = 0, the corresponding states do not feel the perturbation, thus these 
eigenfunctions and the corresponding eigenvalues are identical for the unperturbed and 
perturbed billiards. Next we assume that the perturbed eigenfunctions G obey the 
equation 

{A + k 2 )G(r,R;k) = (2) 

outside of the scatterer of the radius a. In what follows we assume a — > 0. We shall 
see below that G is nothing but the Green function of the unperturbed system. To 
recover an appropriate boundary condition at the perturbation point let us consider 
the asymptotics of the function G(r, R; k) outside of the scatterer when r tends to R. 
Rewriting ([2]) in cylindrical coordinates we obtain 

1 d 2 



1 d d 

pdp^dp p 2 dip 2 



+ k 2 



G(r,R; k) = 0, 



(3) 



where p = |r — R| and <p is the angle between the vector r — R and x-axis going along a 
side of the rectangle. We require G(r, R; k) to be cylindrically symmetric in the vicinity 
of the point r = R, which assumes that the scatterer is cylindrically symmetric. Thus 
we obtain 



1 8 d (2 
p op op 



0. 



(4) 



G(t -»■ R, R; k) 

The solution of the last equation is 

G(r -> R, R; k) = Cl (R; k)J (kp) + c 2 (R; k)Y (kp), (5) 

where Jo and Y are Bessel functions of the first and the second kind respectively and 
C\,C2 are constants parametrically depending on R and k. The last equality should be 
understood in asymptotic sense only, since Jo and Y do not obey the proper conditions 
at the outer boundary of the billiard. 

Using the asymptotic form of Y (z 0) [11] 



Y n (z 



°) = ; 

7T 



+ 7 



0(z 2 \n z) 



where 7 is the Euler constant, and the "identity" 



R|) =5(r-R), 



A(-ln(A;|r 



(6) 



(7) 
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Figure 3. The illustration of the logarithmic singularity of the perturbed eigenfunction 



we find that in the limit a — > the function G(r, R; k) obeys the equation 

(A + A; 2 )G(r,R;A;) = 5(r - R) (8) 

if we assume C2(R; k) = 1/4. Another choice of the constant C2(R; k) would only lead to 
a different normalization. Taking into account the boundary conditions for the function 
G(r, R; k) at the outer boundary of the billiard and (jSJ) we conclude that the perturbed 
eigenf unctions are the Green functions of the unperturbed billiard. 

We are now going to derive the proper boundary condition at the perturbation 
point. First we separate the Green function into its regular and singular part, 
respectively. 

Combining (JSJ), (JSJ) we obtain 



G(t -> R, R; k) = -L In (f) + &(R; k) + 0(z 2 In 



where z = kp, b is some arbitrary length, and 



&(R; fc) = ci(R; fc) + 



2tt 



In 



+ 7 



(9) 



(10) 



Here £&(R; fc) is the renormalized Green function. In figure[3]we illustrate the singularity 
of the Green function near the perturbation point. 

Let us now consider the rectangle with a small pricked circle of radius a whose 
center is situated at the point R. We denote it Q a . Then we consider the linear space of 
functions consisting of two subspaces: (1) the subspace of functions /^'(r, R) vanishing 
at the outer boundary and possessing the asymptotics 



./ 



R) 



B 



2tt \b 



+ 0(z 2 In z) 



[ID 



where B and £ are some constants, and (2) the subspace of regular functions p 2 '{r) 
vanishing at the outer boundary of the billiard such that /^ 2 ^(R) = 0. G(r,H;k) 
belongs to the subspace (1) because of its asymptotic behavior (jHJ). The unperturbed 
eigenfunctions ^(r) belong to the subspace (2). Now we study the action of the the 
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operator —A = 
hermicity gives 



-V 2 on the space of functions defined above. The requirement of 



Iim[</W|-A|/W) a -</«|-A|/W>:] =0, 

a— >0 



(12) 



where /W, f ^' with z, j — 1,2 are arbitrary functions taken from the subspaces (1) and 
respectively, and 

<Prf®*(-Af®). 



(/ W |-A|/W> a 



(13) 



By means of the Green's theorem we obtain 



(/W|-A|/W> B -(/W|-A|/W>: 



2na 



f®*(df®/dp)-f®(df®/dp) m 



(14) 



p=a 



Equation (|T4"j) shows that ffl2|) holds automatically if /W and both belong to the 
subspace (2). Assume now that i = 1 and j = 2 or vice versa. Using the asymptotics (J9]) 
we see that (fT2l again holds automatically for any /W and f( 2 \ Thus the case i,j = l 
implies the only nontrivial condition superimposed by the hermicity requirement of the 
constructed operator. Let us take two functions and f% with asymptotic expansions 



+ 0(z 2 In 2). 



Substituting (JTSj) into flUD we find 



lim 

a->0 



(a (i) i - Ai/n. - (/r i - A|/ro: = sr^r - 6 



0. 



(15) 



(16) 



The equality ([TBI) must hold for any values of £j. This leads to the conclusion that 
for all functions from the subspace (2) the constant £ in (fill is real and t/ie same. Let 
us chose a certain value £ = — D. Then the boundary condition at the perturbation 
point reads: 

£ + D = 0. (17) 

provided that the length b is fixed. Comparing (jUJ) with (ITTj) we find that for the 
perturbed eigenfunctions G(r, R; k) the constant £ is equal to £b(R; fc). Then (TT7T) gives 

&(R; fc) + D = 0. (18) 

Substituting ( flOl) into ( |T8l) we obtain 



ci(R; k) + 



2tt 



In 



~2J 



+ 7 



+ L> = 0. 



(19) 



The proper boundary condition cannot depend on the arbitrary length b, but rather 
should depend on a parameter, characterizing the inner structure of the perturbation. 
Therefore the length b should be canceled in ffl9|) by a proper choice of D. This can be 
achieved by the following choice: 



D —ml? 

2tt 



(20) 



Singular statistics revised 



7 



where (3 is the scattering length of the perturbation. The value of the length (3 can 
not be obtained from the consideration above since we cut out the area containing 
the perturbation, thereby loosing the information on it. Thus we draw the conclusion 
that the scattering length is the only parameter describing the perturbation in the limit 
b C 1. Substituting ( 1201) in (fl8l) and using (flOl) . we find 

0. (21) 



^(R; fc) = Cl (R; k) + -L 

Z7T 



G(r -> R, R; fc n ) = — In f £ ) + 0(;z 2 In 2). (22) 



The perturbed part of the spectrum may now be obtained from the solutions k 2 = k 2 n 

of flrij. 

Replacing b by /3, k by fc n in (jUJ) and using the equality £/?(R, k n ) = we find the 
asymptotic expansion of the perturbed eigenfunction corresponding to the eigenvalue 

d(r -y R R k-\ = 

2vr \/3, 

The leading term of the asymptotics (1221) becomes zero when p = (3. This fact can be 
used to determine experimentally a scattering length of a given perturbation. 

We note that the proper definition of the scattering length was missing in Refs. 
[HIE]. This has led to the deficiencies discussed above. However in the monograph [8], 
devoted to point perturbations, the scattering length in two-dimensional problems was 
properly introduced. 

Several conclusions on the level-spacings distribution can be drawn already from 
(l2Tj) . Indeed, from (jHJ) and (ITU]) we conclude that ci(R;/c) has the same poles as the 
Green function of the unperturbed billiard. From (12T1) we obtain that when k tends 
to infinity, the eigenvalues of the perturbed billiard approach the eigenvalues of the 
unperturbed one. Indeed, close to the eigenvalue k n the function ci(R, k) may be 
approximated by const/(/c 2 — k^) whence follows: 

f ' OUSt + V ,(k)=0, #)4 h f?)» 1 - ( 23 ) 



k 2 -k 2 n ,pv ' ' ,pv ' 2tt V 2 , 

Then k 2 — k 2 = —2n const/ ln(fc/3/2). In the limit k — » 00 we recover the original 
spectrum of the billiard! Since the statistics of inter-level spacings for the unperturbed 
rectangular billiard with chosen side ratio is Poissonian [2], we conclude that the same 
statistics for high-lying eigenvalues of the Seba billiard is also Poissonian. This is the 
most important conclusion of the paper, which contradicts the prediction given in [JJ. 



3. Ewald's representation of the renormalized Green function 

For an explicit calculation of the spectrum of the perturbed billiard from (j2T!) an 
expression for ^(R; k) is needed, which for the general cases by no means is a trivial 
task. For the rectangle it can be obtained by an application of Ewald's method. The 
derivation is technical, and anybody not interested in the details may proceed directly 
to 
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Figure 4. The illustration of the sources distribution used in the image representation 
of the Green function 



In the paper [7] in (|2ip the logarithmic dependence of the dimensionless scattering 
strength T]p{k) was missed. Therefore the proper spectral statistics of the Seba billiard 
still has to be computed. In this section we explain the numerical procedure, based on 
Ewald's method (131 [HI [151 EH CEO US], which allows us to compute the renormalized 
Green function and then, from ([21]) . the spectral statistics. 

We start from the eigenfunction representation of the Green function for the free 
billiard 

inn 

k 2 -E n 



G(r, R; k) = ^ 2^ h 2 _ R ' ( 24 ) 

where R = (x', y' 



n=l m=l 



^ m{x ' v) = ^ sm IT ) sm ^lf) ' (25) 



irnx 



\ 2 / \ 2 

irn\ ( Tim \ 



E ™ - UJ + UJ • (26) 

<i x and d y are the two sides of the rectangle. When x — > x' and y — > y' then the series (J24l) 
diverges logarithmically. This is just another manifestation of the well-known singularity 
of the Green function for r —>■ R, see ([9]), which is a local feature and does not depend 
on outer boundary conditions. This suggests that the eigenfunctions representation is 
not the appropriate choice to compute the renormalized Green function, but that the 
images representation [16] might be preferable. 

The images representation of the Green function reads 

oo 1 

G(r, R;k) = E (-l) Sl+S2 C/(r, R SlS2 + R nm ; k), (27) 

n,m=— oo si,S2=0 

where 



R Sl s 2 = [{-l) Sl x', (-I) 52 /) , R« m = (2nd x , 2md y ), 



(28) 



and Gf(r,H;k) is the Green function for the two-dimensional plane (see figure H]). 
The images representation is much better suited to compute the renormalized Green 
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function, since when r tends to R, the term n = m = 0, S\ = S2 = is the only one 
in the representation logarithmically tending to infinity. Now the divergency can be 
subtracted analytically However the image representation does not solve the problem 
yet since it converges absolutely only if Irak > 0. To overcome this obstacle Ewald 
[T3] proposed the dual representation keeping features of the images as well as the 
eigenmodes representation. 

Below we follow the works [HI [15], [16] . Let us first find a convenient representation 
for Gf. To this end we consider the following initial- value problem for the function 
g(t;r,R, k) 

^ = (A + k 2 )g, g(t = 0;r,R,k) = 5(r-R). (29) 
Then Gf can be written as 

G f (r,R;k) = - [ g(t;r,R,k). (30) 
Jc 

The contour C should start at t — and tend to infinity in such a way that g tends to 
zero. Obviously 

g(t;r,R,k) = e k2t K(t ]r ,R), (31) 
where the heat kernel K(t; r, R) can be found by a separation of variables 

K(t ] r,R) = —e~^ 2 ^ t l (32) 
Airt 

Finally we obtain 

^•^-IIM^-^)- (33) 

The simplest contour of the integration is the imaginary half-line going from zero to ioo 
(contour C\ in figureEl(a) ). Using this contour we conclude ([IT]. Entry B. 187(2)) that 

G / (r,R;A;) = -^ 1) (A;|r-R|), (34) 

where Hq is the Hankel function of the first kind. However, the asymptotic behavior 
becomes hidden in this representation. To recover the asymptotic behavior we will 
use another contour of the integration. We see that for small values of \t\ the best 
convergency provides an interval lying on the real axis from zero to some positive value 
tEw (see figure [5] (a) ). We will call this value the Ewald parameter. From the other 
side the best convergency for large values of \t\ would provide the half-line going from 
some negative value (we choose it to be equal to —tEw) to —00. What remains is to 
connect these parts to make a contour. We connect them by a half-circle C^- Finally 
the constructed contour is equivalent to G\ since the integrals along the quarter-circles 
C 3 , C 4 (plotted by dashed lines in figure [5] (a) ) tend to zero when the radius of C 3 tends 
to infinity and the radius of C 4 tends to zero. 
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Figure 5. The contours of the integration for the free Green function Gf (a) and the 
contour of the integration for the Ewald's representation of Gf (b) 

Using the constructed contour we can obtain the asymptotics of the type ([9]) from 
the representation ( 1331) . The developed technic will be used further to compute the 
renormalized Green function. We write 

(35) 



where 



G f (r, R; k) = Gf\r, R, k) + Gf\r, R; k), 



Gf\r,R; k) 



) 

g(t; r, R, k)dtj 



P P — oo 

Gf\r,R;k) = - g{t;r,R,k)dt- g(t;v,R, 

J C*2 J —tEw 



k)dt. 



(36) 



Introducing the notations u = (r — R) 2 /(AtEw), v — k 2 tE W , we write G^ as follows 



G^(r,R; k) 



1,1 dt 



vt—u/t 



(37) 

Jo Ant 1 ; 

When r tends to R, then u tends to zero. To compute the asymptotics of (!37j) for u — > 
we make the following transformations: 



Gf(r - R, R; k) = - jT ^ (e* - 1 + l) e"^ 



4vrt 



~ / ™ (i _ e «) _ 



r — r i - ^ - z 1 ie- - r — 

rl dt 
Jo 4^rt 



/>oo 



~ / J!L(i_ e -t)_ 



Ant 

dt 
l 1 ~Axt 



! 4vrf 

*+ T— 

7o 4vrt 



™ (l_ e *) + J_ lnw 



= ^ + 7 + In u] 



where 



9o(v) = T T (1 " **) , 7 = r 

Jo r ■/ o 



(38) 
(39) 



ml dt„ _ u f°°dt _ t 

(see [TT] (Entry 8.367.12)). Since go(^) has no singularity at t — it can be written as 

- 1 dt 



9o(v) = f j (1 - e-) - / f (1 - e") 

J */ Co 



(40) 
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where Cq is a half-circle of the unit radius. 



(2) 

Since G y has no singularity when r — > R, the leading term of its asymptotics is 

Gf(W = -/ * -_^-* e -. (41) 



Using (EHD, (30}, (SI1) we find 

^ 1} (r -> R, R; jfe) + (R, R; jfe) 

The last equality does not depend on the choice oft Ew . To prove it we can take v as the 
independent parameter, then u = k 2 (r — R) 2 /(4t>). Differentiation of the last equality 
with respect to v gives zero. The reasonable choice of v should not lead to exponentially 
large values of go(v). Indeed due to fj4"2]) such a large contribution is somehow artificial 
since it is annihilated by Gj . Thus it is natural to take v — 1. Then (14*21) gives 



G f (r^R,R; k) ~ i- 



, . &|r - R. . 
ln i 2 1+7 



This calculation has demonstrated that it is important to divide the free Green function 
in two parts: Gp and G^J\ The first part describes the space singularity, and the 
second part makes the contribution into the regular part of the asymptotics. 

Let us turn now to the image representation (1271) of the Green function of the 
rectangular cavity. We again write the free Green function Gf in the form 

G f (r,R;k) = Gf\r, R; jfe) + Gf \r, R; jfe), (44) 
where Gf is defined by (13T)) and 

G ?^=-L£M et - {1 ^)- (45) 

The best choice of the contour C5 will be discussed later. We shall see that the choice 
used to compute the asymptotics (l4"3l) does not fit. Substituting f[44|) into (12T1) we obtain 

G{t, R; jfe) = G (1) (r, R; jfe) + G (2) (r, R; jfe), (46) 



where 



G«(r,R;fc)= Yl E (-l) n+8 'Gf\T,R ai9a +R nm ',k), (47) 

n,m=— 00 si,S2=0 
00 1 

G( 2 Hr,R;fc)= ^ ^ (-l^Gf (r, R S1S2 + R nm ; k). (48) 



n,m=— cxd si,S2=0 

To improve the convergency of series for G^ we use the identity 
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which can be proved by applying the Poisson sum rule [H] to the function 

f(x) = e-*K* t \ (50) 
Performing the resummation we obtain 
G (2) (r,R; k) 

oo 1 „ 

'y ^ ^iimx / d x +i-Kmy / d y 'y ^ ^_^yi+S2 J e [k 2 -(nn/d x ) 2 -(TTm/d y ) 2 ]t^ (51) 



4d x dy j Q 

y n,m=—oo si,S2=U s 

Now we see that the integral over C5 should converge for positive as well as for negative 
values of the real part of k 2 — (Trn/d x ) 2 — (nm/dy) 2 provided that Imk > 0. Therefore 
we have to assume t — > zoo along the contour C5. This leads to the choice of the contour 
shown in figure [5] (b). 

Performing the integration we obtain 

1 ~ [k 2 -(-Kn/d x ) 2 ~^ m /d y ) 2 }t Ew 

^d x d y n k 2 - (nn/d x ) 2 - (nm/dy) 2 

1 

x ^_^y 1 +S2 e -iirn(-l) s ix' /d x -iirm(-l) s 2y' /d y ^g2) 

S1,S2=0 

Summarizing over si, S2 we get 

J2 ^ly^^ni-ir^'/d^mi-ir^/dy = _ 4sin sin /'^lA (53) 

Sl,S2=0 ' \ » / 

Now we can perform summations over n and m: 
1 



G (2) (r,R; fc) 



d x dy 

!Ll^ Urn I — " I Jimx/dx+mmy/dy 



E( imx' \ ( irmy' \ 



x > sm I — ; — I sm I — ; — I — ; t^t^ ; rrr^e 

hi 



n,m= — 00 



k 2 — (im/d x ) 2 — (irm/dy^ 2 



E 



4 v / 7ma;\ / 7rmy\ ( imx \ ( ixmy 
sin — - — sin — - — sin — - — sin 



d x dy \ d x ) V dy J V d x J V 

a n,m,=l x ' \ a / x / \ a 

e [fc 2 -(7rn/^) 2 -(7rm/d y ) 2 ]t Btl , 

A; 2 — (nn/d x ) 2 — {irm/d y ) 2 ^ 

Formulas (1371) . (1471) and ( 1541) give the Ewald representation of the Green function for the 
rectangular billiard. The integral in (l37j) has to be computed numerically. Now we can 
recapitulate the advantages of the Ewald's method. First of all, both series and 
are exponentially convergent. Thus we can take the analytic continuation and choose 
real k. Second, we have separated the part G^ responsible for the space singularity 
from the part G^ responsible for the poles information. Indeed, G^ exponentially 
converges even when r = R. Third, only in G^ 1 ' there is the term corresponding to 
Si = S2 = 0, n = m = which asymptotically tends to infinity when r — > R. The rest 
of the series is exponentially convergent. The last observation allows to compute the 
renormalized Green function. 
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Though the Poisson resummation is a common tool used to get the Ewald's 
representation of the Green function [T3l [Til [T51 [T6l [18], one can avoid it and obtain 
formulas ( 1471) and ( |54|) easier (see the Appendix for details). 

Using the asymptotic expansion (138)) we obtain the ea;ac£ Ewald representation for 
the renormalized Green function: 

6g(R; fc) = ^[3o(fcV«)+7] + 

oo 1 

+ X] (l-^,o5 m ,o^ 1 A 2 ,o)(-l) Sl+S2 ^ 1) (R,R SlS2 +R nm ;A;) + 

n,m=— oo si,S2=0 

+ i ln (il) +G(2,(R - Ri ' : »- < 55 » 

with G^p to be computed numerically from the integral ( 1371) . and G^ from the sum 
(151)1 . Now we can compute the perturbed part of the spectrum from the condition 
£/3(R, fc n ) = 0, see (ED), using the representation ( 1551) . This final equation in contrast 
to (1T9l) lost the clearness since it depends on the as yet not defined Ewald parameter 
tEw To define it we first consider large values of k. Then to avoid exponentially large 
values of the function go(k 2 tE W ) as well as exponentially large amplitudes of terms with 
small numbers n, m in the expansion of G^ we put tE W — l/k 2 . Obviously this choice 
is inappropriate for k — > 0, since this would mean to compute a huge number of terms 
in So, finally the Ewald parameter can be chosen as 

/ 1/k 2 , if k> ko, . , 

tEw ~ \ l/kl if k < fc , (56) 

where k\ = (ir/d x ) 2 + (n/dy) 2 is the lowest eigenvalue of the unperturbed system. To 
investigate the spectral statistics we can assume tEw = 1/& 2 . Then 02 II) reads 

^>(l)+7] + 

oo 1 

+ Yl (l-^, 5 mi 0^ 1 A 2 ,0)(-l) Sl+S2 4 1) ( R ' R ^2+Rn m ;fc) + 

n,m=— oo si,S2=0 

+ 27 ln (^) +G(2)(R ' R;A;) = - (57) 

Now (1571) resembles (TT9|) . so the main conclusions made above could be repeated. The 
equation (1571) is alike (3) in [7J , apart from the fact that G^ 2 ' (R, R; k) is not a finite sum 
and the rest in (157)) is not a polynomial as a function of k 2 . Equation (1571) is exact and 
especially fits for the numerical study, since it contains exponentially convergent series. 

For large k the double sum in (157)) can be neglected and the rest looks very similar 
to the "TV-poles" approximation [10] . Indeed in this regime the spectrum of the billiard 
can be found from the equation 

'k(3\ , 7 + <7o(l) 



G( 2 )(R,R;fc) + ^[ln (M) + 



0, (5* 
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Figure 6. Graphical interpretation of (|57l) . Vertical grey dashed lines correspond to 
eigenwavenumbers of the unperturbed billiard. Thin solid line shows (R, R; k) as 
a function of k. Thick solid line corresponds to the remainder in (|57p taken with minus 
sign. In the figure = 1. 

where only a finite number of terms in the expansion of can be taken into account 
due to the exponential convergency. Figure [6] shows a graphical interpretation of fl57|) . 
In our calculations we found that approximation ( 1581) works perfectly above the first 
resonance already. 

4. Integrated density of states 

In what follows we are interested in level-spacing statistics for the subset of perturbed 
eigenvalues of the Seba billiard. There are several reasons to restrict ourselves to 
the statistics of the subspectrum. First of all the influence of the perturbation is 
more pronounced if one considers only the perturbed part of the spectrum. This is 
probably the reason why in the pioneering work [1] only the statistics of the subspectrum 
is considered. Another reason to consider subspectrum's statistics is ( |2TT) . which 
determines only the perturbed subspectrum. The graphical interpretation of fl57|) gives 
already an idea on the structure of the perturbed subspectrum (see figure [6]), while 
considering the unperturbed subspectrum as well we loose the clearness. The last reason 
to consider the statistics of the perturbed subspectrum only is the direct correspondence 
of the perturbed subspectrum to the spectrum obtained from the reflection measurement 
with a single antenna introduced at the point of the perturbation. In such an experiment 
the unperturbed subspectrum is not seen at all since the corresponding eigenstates, 
vanishing at the perturbation point, can not be excited. 

If one considers only the perturbed subspectrum it makes a difference whether the 
ratios x'/d x and y'/d y are rational or irrational numbers (see figure [7]). The difference 
arises from the fact that for irrational numbers all eigenfunctions are perturbed while 
for rational ones a part of eigenfunctions remains unperturbed. 

To compute the statistics and compare it with GOE, semi-Poissonian and 
Poissonian predictions one should first unfold the spectrum to a mean level spacing 
of one. This can be achieved by the following definition of the scaled eigenvalues: 

= NQti, (59) 
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-dx/qi 



dx 



r " ! 



1 



dy 



• <it pi/qi 



Figure 7. "Subbilliards" corresponding to rational ratios x'/d x , y'/d y . The point 
(x', y') is shown by the black disk. 

where N(z) is a smoothed function counting a total number of eigenvalues k\ less 
then z, i.e. the integrated density of states. If a spectrum of a system is known, the 
function N(z) can be obtained from a numerical fit. For a conventional unperturbed 
two-dimensional billiard one can use the Weyl estimation of the integrated density of 
states (see e. g. [19J) 

N w{z) = ^ z ~ + A w , (60) 

where A\ is the area of the billiard, A 2 is its circumference and Ay/ is a constant. 
For the unperturbed rectangular billiard with the sides d x , d y we obtain A\ = d x d y 
and A2 = 2(d x + dy). While the Weyl estimation holds for the whole spectrum of the 
unperturbed billiard it can not be directly applied to its subspectrum as well as to the 
perturbed subspectrum of the Seba billiard. However to fit numerically the integrated 
density of states one can still assume that the function to be found has the form (1601) 
with some unknown coefficients A\, A 2 , Aw- 

The scaled level spacing corresponding to the nearest eigenvalues k\ and k^ +1 is 

s n = - = N(k 2 n+1 ) - N{kl). (61) 

Thus the constant term Aw in ( IBTJl) does not influence the statistics. The mean level 
spacing 

1 - 1 

< s ) = M E s « = M [N{kli+l) ~ N{kl)] ~* 1 (62) 

71=1 

when M — > 00 in accordance with the rescaling requirement. 

Though the integrated density of states corresponding to the perturbed 
subspectrum can be fitted numerically, it is possible to estimate it a priori. Indeed in 
figure M one sees that between two successive eigenvalues of the unperturbed spectrum 
corresponding to poles of the function there always exists an eigenvalue of the 
perturbed billiard. Thus the number of perturbed eigenvalues of the billiard below 
z should coincide (up to a single eigenvalue) with a number of eigenvalues of the 
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unperturbed billiard below z corresponding to nonvanishing eigenf unctions at the point 
of the perturbation. Obviously this result does not depend on the value of the scattering 
length. 

Following the argument given above we can compute the expected function N e (z) 
right from the unperturbed billiard, where N e (z) is equal to the integrated density 
of those states whose eigenfunctions do not vanish at the perturbation point (x',y'). 
Let us assume that x' = d x pi/qi, y' = d y p 2 /q2, where p±/qi and P2/92 are irreducible 
fractions. From figure [7| one can draw the conclusion that eigenfunctions of the small 
hatched "subbilliards" with Dirichlet conditions at all boundaries are eigenfunctions of 
the initial billiard and vanish at the point (x',y f ). According to the Weyl formula the 
number of eigenvalues below z can be estimated as 

N v (z) = ^z - 4/9 1 + ^ + A v , (63) 

4lTq 1 Z7T 

N h (z) = j^z - dx \ dy,q2 ^-z + A h (64) 

for the vertical and horizontal hatched billiards respectively. Here A v and Ah are 
some constants. We have computed twice the eigenvalues of the billiard obtained as 
an intersection of these subbilliards. Its number of eigenvalues can be estimated as 

at 1 \ d x d y d x /qi + d y /q 2 r , 

N vh {z) = - z \/z + A vh . (65) 

47rgig 2 27r 

Finally, the number of eigenvalues corresponding to vanishing eigenfunctions is 
N v (z) + N h (z)-N vh (z) 

- + --—)^z-^±^yG + A v + A h -A vh . (66) 
Subtracting the last estimation from the total number of eigenvalues below z 

N w (z) = ^l z - ^±A^ Z + Aw (67) 

we obtain the following estimation for the number of eigenvalues corresponding to 
nonvanishing eigenfunctions: 

N e (z) = ( 1 - - - - + d -^z + A e , (68) 

where A e = Aw + A v h — A v — Ah- Surprisingly the surface contribution ~ sfz vanishes. 

In figure M the estimation (1681) and the numerical fit of the form (1601) for the 
integrated density of states are shown for the comparison. One can see that the 
estimation fl68|) works very well. 

From ( 1581) . ( 1681) one can compute a shift of a resonance induced by a point 
perturbation as compared to the mean level spacing. Indeed when k — > E nm the function 
G^(K, R; k) (El) tends to the following expression: 

G^(K,K;k) _ 4 sin 2 ^)^^) (g9) 
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Figure 8. Integrated density of states of the perturbed subspectrum of the Seba 
billiard. Figure (b) represents the zoom of figure (a). The perturbation with the 
scattering length j3 = 1 is placed in the center of the billiard. Stepwise line corresponds 
to the computed subspectrum, dashed line corresponds to (|68|) . smooth solid line 
corresponds to the numerical fit of the form (|60|) . 



From (IB"8l) we find the mean level-spacing (AE): 

(AE) = 1/N'(z) = -j-j-, 1 = 1-1-1 + J-. (70) 
d x d y Q qi q 2 qiq 2 

Substituting ( l69l) in ( J58l) and using (|70l we find for the relative shift of the resonance: 
k 2 - E nm _ d x d y 2 4sin 2 (7rna: / /4)sin 2 (7rmy7^) 
(AE) ~ AttQ [ nm) = Q[\n(E nm (3yA) + 7 + <fo(l)] ' 1 ' 

Let us estimate the number of resonances needed to show the transition to the Poissonian 
level-spacing statistics. Then the relative shift should be very small for all sufficiently 
large numbers n and m. The sufficient condition is 

|[ln(£ nm /3 2 /4) + 7 + g (l)} » 1. (72) 

Depending on the values of Q and (3 the value of E nm can be very large. 



5. Level-spacing statistics 

In this section we are presenting a number of numerical results. Figures [9lfTT1 show level- 
spacings distributions for the perturbed part of the spectrum for various situations to 
be discussed in detail below. For comparison the curves corresponding to Poissonian, 
semi-Poissonian and GOE distributions are plotted by solid, dashed, and solid lines, 
respectively. In each of the figures subfigure (a) shows the level-spacings distribution 
for the unperturbed system to make sure that the distribution is really Poissonian, 
since it is well-known that there may be deviations for small distances depending on 
the side ratio of the rectangle. Subfigures (b)-(d) show level-spacings distributions for 
(3 = 1, 0.1, 0.02 respectively. Note that a decreasing value of (3 means an increase of the 
perturbation. 

In figure [9] the scatterer is placed in the center whereas in figure [10] it is at the 
point (0.55^,0.65^). In both cases about 1500 lowest perturbed eigenvalues have 
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Figure 9. The level spacings statistics for the first 1617 resonances. The subspectrum 
of the unperturbed billiard corresponding to nonvanishing eigenfunctions at the center 
of the billiard (a). Perturbed subspectrum of the Seba billiard with a perturbation 
placed at the center, (3=1 (b). (c) the same as in (b) with (3 = 0.1. (d) the same as 
in (b) with = .02. 



been considered. For (3 = 1, i.e. for a weak perturbation, the distribution shows a linear 
repulsion for small distances and an exponential tail, but not a semi-Poissonian behavior 
in the strict sense (dashed line). With (3 = 0.1, 0.02 there is a gradual transition to a 
broader distribution, resembling GOE one for the scatterer in the center (figure [9](d) ) 
and a semi-Poissonian distribution for the scatterer in the off-center position (figure 
ITOTd) ) . This observation would deserve more quantitative treatment, but this goes 
beyond the scope of this paper. Qualitatively it may be understood from the fact that 
in the center of the billiard all perturbed eigenfunctions have the same value, while for 
an off-center position there is a distribution of eigenfunctions amplitudes giving rise to 
a corresponding distribution of resonances shifts. 

Figure [11] finally shows the level-spacings distribution again with the perturbation 
in the center but now for the numbers of perturbed eigenvalues from 25000 to 27000. 
Comparison of figures [9] and [11] shows a pronounced change of the distribution towards 
Poissonian with increasing eigenvalues numbers. This is particularly evident for the 
weaker perturbations (3 = 1, 0.1 (figure fTT1(b)-(c) and demonstrates the main result of 
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(c) (d) 

Figure 10. The level spacings statistics for the first 1482 resonances. The subspectrum 
of the unperturbed billiard corresponding to nonvanishing eigenfunctions at the 
position (0.55g? x , 0.65^) of the billiard (a). Perturbed subspectrum of the Seba billiard 
with a perturbation placed at the point (0.55(4;, 0.65d y ), [3 = 1 (b). (c) the same as in 
(b) with P = 0.1. (d) the same as in (b) with j3 = .02. 



this paper: with increasing eigenvalues numbers eventually the lev el- spacings distribution 
of the unperturbed system is recovered. 

6. Conclusions 



Let us now conclude. First of all we have presented in the paper the complete solution 
of the spectral problem for the rectangular billiard with a single point perturbation. In 
contrast to previous studies [H [7j we have shown that the statistics of the Seba billiard 
tends to a Poissonian when the number of levels taken into account tends to infinity. 
The estimation given at the end of Section 5 showed, however, that the transition to 
Poissonian statistics appears, depending on the scattering length, only at exponentially 
large quantum numbers. The solution is based on the Ewald representation of the 
renormalized Green function ( 1551) . This representation contains exponentially rapidly 
convergent series. Together with the Ewald representation of the usual Green function 
(157)1 . (1471) . (1541) the presented approach is a powerful tool to analyze various experiments 
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(d) 



Figure 11. The level spacings statistics for the resonances 25000-27000. The 
subspectrum of the unperturbed billiard corresponding to nonvanishing eigenfunctions 
at the center of the billiard (a). Perturbed subspectrum of the Seba billiard with a 
perturbation placed at the center, (3=1 (b). (c) the same as in (b) with /3 = 0.1. (d) 
the same as in (b) with j3 = .02. 



made in rectangular billiards. 
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Appendix 

In the Appendix we show that the explicit application of the Poisson resummation in 
the derivation of Ewald's representation of the Green function can be avoided. These 
findings simplify the technical calculations. 
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Let us consider the initial problem (I29p in the rectangular billiard with proper 
boundary conditions. Then the solution can be written in two equivalent forms: in 
the form of the images representation g>j(t; r, R, k) and in the form of eigenmodes 
representation g e (t; r, R, k): 

oo 1 

9l (t;r,R,k) = Yl (-^) sl+S2 9{t;r,R slS2 +R nm ,k) (A.l) 



n,m=— oo si,S2=0 
oo 



g e (t;r,R,k) = -^- e^-^^-^^V-m^y^mn^.l/)- ( A -2) 

a * a v n ,m=l 

Then we write 

G(r, R; k) = G (1) (r, R; k) + G (2) (r, R; k), (A.3) 

where 

(r, R; fc) = - / gi (t; r, R, k)dt, (A.4) 

G (2) (r, R; k) = - [ g e (t; r, R, k)dt. (A.5) 

Performing the integration in Eqs. ( 1A.4I) . f lA.5[) we immediately get ( HTj) . ( |54l) . 
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